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ABSTRACT 

Context. The energy released during solar flares is believed to be stored in non-potential magnetic fields associated 
with electric currents flowing in the corona. While no measurements of coronal electric currents are presently available, 
maps of photospheric electric currents can now be derived from SDO/HMI observations. Photospheric electric currents 
have been shown to be the tracers of the coronal electric currents. Particle acceleration can result from electric fields 
associated with coronal electric currents. We revisit here some aspects of the relationship between particle acceleration 
in solar flares and electric currents in the active region. 

Aims. We study the relation between the energetic electron interaction sites in the solar atmosphere, and the magnitudes 
and changes of vertical electric current densities measured at the photospheric level, during the X2.2 flare on February 
15 2011 in AR NOAA 11158. 

Methods. X-ray images from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) are overlaid on 
magnetic field and electric current density maps calculated from the spectropolarimetric measurements of the Helio- 
scismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO) using the UNNOFIT inversion and 
Metcalf disambiguation codes. X-ray images are also compared with Extreme Ultraviolet (EUV) images from the SDO 
Atmospheric Imaging Assembly (AIA) to complement the flare analysis. 

Results. Part of the elongated X-ray emissions from both thermal and non-thermal electrons overlay the elongated 
narrow current ribbons observed at the photospheric level. A new X-ray source at 50-100 keV (produced by non- 
thermal electrons) is observed in the course of the flare and is cospatial with a region in which new vertical photospheric 
currents appeared during the same period (increase of 15%). These observational results are discussed in the context 
of the scenarios in which magnetic reconnection (and subsequent plasma heating and particle acceleration) occurs at 
current-carrying layers in the corona. 

Key words. Sun: flares - Sun: particle emission - Sun: X-rays - Sun: magnetic fields - Acceleration of particles - RHESSI 
- SDO/AIA - SDO/HMI 


1. Introduction 

It is now commonly admitted that solar flares are the re¬ 
sult of the sudden release of magnetic free energy in the 
corona, this free energy being stored in non-potential mag¬ 
netic fields associated with electric currents flowing in the 
corona. As there are currently no measurements of the vec¬ 
tor magnetic fields in the corona, there are also no measure¬ 
ments of coronal electric currents. Electric currents have 
been only determined in a few cases at the photospheric 
level using vector magnetic field measurements achieved 
in photospheric lines. Whi le some pioneering work was 
done in the 1970’s by e.g. iMoreton fc SeveriM (| 19681) at 
the Crimean Astrophysical Observatory, many observations 
of vector magnetic fields and related electric currents have 
been obtained later in the 1980’s and 199Q’s with ground- 
based vector magne tographs (see e.g. iHagvard et al.|[T9^ 
ICanfield et al.|[T9^ . Nowadays, the polarimetric measure¬ 
ments obtained continuously, and with an unprecedented 
cadence and spatial resolution with the Helioseismic and 
Magnetic Imager (HMI) on Solar Dynamics Observatory 
(SDO) allows to obtain maps of photospheric vertical elec¬ 
tric currents for flaring active regions and even more to 


study the evolution of these currents on a time resolution o f 
12 minutes (see e.g. lPetri3l2Q12L [^QljUJanvier et al.ll2Ql4 ). 

The relation between the energetic phenomena (plasma 
heating and particle acceleration) occurring in solar flares 
and the electric current systems in active regions has been 
the subject of many studies for several decades. In particu¬ 
lar, the link between the photospheric electric currents and 
the energetic electron precipitation sites has been investi¬ 
gated in several events, originally to discuss the relevance 
of the different flare models presented in the literature. The 
first studies used Ha observations to characterize electron 
precipitation sites since no HXR imaging observations were 
available. Using data from the Crime an Astrophysical Ob¬ 
servatory, iMoreton fc Severn^ (|1968[ ) examined the spatial 
relationship between Ha flare kernels and maxima of verti¬ 
cal electric current densities. Eor 80% of the 25 events stud¬ 
ied, a spatial coincidence (< 6”) was found between the loca¬ 
tion of the center of the bright Ha kernel and the location of 
the strongest electric current densities (> 8 mA/m^). Such 
a spatial coincidence between Ha flare kernels and strong 
elec tric current densiti e s was later confirmed by other stud¬ 
ies. iLin fc Gaizanska!^ (|1987[ ) found in particular that the 
sites of strongest Ha emissions were cospatial with the ver- 
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tical current systems at the photospheric level, to within 
2”. However, it was also found that even if some of the flare 
kernels were close to the electric c urrents, many of them ap - 


sar 


^1199011 


Leka et H] 


peared near the edges o f currents (iRonianov fc 
In a s eries of papers, ICanfield et al.l ( 19931) 

(|1993[ ) ; Ide La Beauiardiere et al.l (|1993[ ) revisited the ques¬ 
tion of the relation between electric current systems and 
electron precipitation sites using vector magnetograph data 
from the Mees Solar Observatory and information from Ha 
line profiles which allow to disentangle Ha signal from elec¬ 
tron preci pitation and from high pressure. Studies of sev- 
eral flares (|Leka et al .11 19931: ide La Beauiardiere et al.lll993[ ) 
confirmed that sites of intense nonthermal electron precip¬ 
itation do not coincide with the regions of strongest verti¬ 
cal currents at the photosphere, but tend to occur on the 
shoulders of channels of high vertical current density rather 
than at the vertical current density maxima. In a further 
study, iDemoulin et al.l (|l997l ) found for several events that 
flare ribbons were located in the vicinity of strong elec¬ 
tric currents in the photosphere. They furthermore derived 
the locations of the photospheric footprints of regions of 
rapid change in magnetic line connectivity called quasi- 
separatrix layers (QS Ls). They found as in other observa¬ 
tiona l studies (see e.g [Machado et al.|[T9^ iMandrini et al.l 
1 19951: iBagala et al.l Il995f ) that flares tend to occur near 
the locations of these QSLs. These QSLs are now known 
to be regions where strong current densi ties can develop 
and where reconnection can occur (see e.g. IDemoulin et al.l 
1 19961: lAulanier et al.l l2QQ5l ) thus establishing a clear link 
between flare energy release, electric currents and magnetic 
reconnection. 


With the arrival of HXR imagers, the link between 
vertical electric currents at the photosphere and electron 
precipitation sites has been revisite d using in pyticu- 
lar YOHKOH /HXT observations f see I Canfield et anil992l: 
iLi et al.lll997[ ). The latter study based on six events con¬ 
firms that HXR emission and thus electron precipitation 
sites are not exactly co-spatial with regions of highest ver¬ 
tical current densities at photospheric levels but are rather 
adjacent to the curr ent channels, therefore confirmin g the 
results obtained by Ide La Beauiardiere et al.l (Il993l) with 
Ha observations. In a later paper. lAschwanden et al.l (| 19991) 
showed how this offset between the maximum of vertical 
currents and the HXR loop footpoints could be explained 
in the context of the 3D geometry of quadrupolar recon¬ 
nection. It must finally be noticed that in the earlier stud¬ 
ies, vector magnetograms (and thus electric current maps) 
were sparsely derived and that the integration time to de¬ 
rive such information could be more than one hour. Several 
hours may also separate the time of the flare and the time 
of the magnetic field measurements. Nowadays, the combi¬ 
nation of polarimetric measurements continuously obtained 
with the Helioseismic and Magnetic Imager (HMI) on Solar 
Dynamics Observatory (SDO) and of HXR observations of 
solar flares obtained with the Reuven Ramaty High Energy 
Solar Spectroscopic Imager (RHESSI) allows to compare 
high quality vector magnetic field (and electric current) 
maps at a cadence of 12 minutes and HXR observations 
of solar flares at exactly the same time. This provides new 
possibilities to study more systematically and in more de¬ 
tails the relation between the location of HXR sources (and 
electron precipitation sites) and the vertical electric current 
densities. 


^ . In clo se relationship with 
l| (|2Q14l) quantified the ver- 


This paper presents the first results of such a study for 
the flare of February 15, 2011. The simultaneous evolutions 
of electric currents and of HXR emission sources during 
the flare are furthermore examined for the first time. The 
February 15, 2011 event is quite appropriate for such a first 
study. Indeed, this major flare (GOES class X2.2) occurred 
in the active region AR11158 when the active region was 
near disk-center. The derivation of vector magnetic field 
maps and of vertical electric current densities requires in¬ 
deed observations of an active region located near the center 
of the solar disk to have good measurements of the vector 
magnetic field at the photospheric level and to avoid projec¬ 
tion effects. The February 15, 2011 flare which is the first 
X-class flare of cycle 24 is the subject of num e rous stud- 
ies since the original analysis by ISchriiver et al.l (12011 ) (see 
e.g. for some of the most recent papers llnoue et al.l 20l4 
for MHD simulatio ns of the flare) 
the present paper, I Janvier et al.l 
tical currents at the photospheric level and studied their 
temporal evolution during the flare. They concluded in par¬ 
ticular that the evolution of both current and flare ribbons 
is related to the evolution of the magnetic field during the 
flare in the context of a 3D standard flare model. The ac¬ 
tive region itself and i ts activity has be e n also extensively 
studie d . In particular. ISun et al.l (l2Q12|j: IVemareddy et al.l 
(|2Q12l ): lAschwanden et al.l (|2Q14l ): IZhao et al.l (|2Q14 j real¬ 
ized different estimations of the e.g. free magnetic energy 
and dissipated energy for several events occurring in this 
active region. 

In section [2] of this paper, we present briefly the dif¬ 
ferent instruments used in this study as well as the tech¬ 
niques necessary to obtain maps of vector magnetic fields 
and electric currents from the measurements. In section [3l 
X-ray and EUV observations of the flare are compared with 
the vertical electric current densities. Section 0] discusses 
the observational results and presents our interpretation. 
A summary of the paper and of the conclusions is finally 
presented in section [5l 


2. Magnetic field measurements and Flare 
observations 

2.1. Magnetic Field Measurements and Electric Current 
Calculations 

The magnetic field and electric current density maps are 
derivated from the data of the Helioseismic a nd Magnetic 
Imag er HMI on SDO (|Scherrer et al.l l2Q12l: ISchou et al.l 
12012 ). The HMI instrument provides images of the entire 
Sun in six narrow spectral bands in a single iron line (Fel 
617.33 nm) and in four different states of polarization. This 
set of 24 images provides spectropolarimetric data to en¬ 
able the calculation of the full vector magnetic field (in¬ 
verse problem), at the altitude of formation of the spectral 
line, which is at the photospheric level in this case. The 
three components of the photospheric vector magnetic field 
can be calculated together with the vertical component of 
the electric current density (see e.g. figure [1] bottom and 
iJanvier et al.ll2Ql3 for further details). However, because 
the spectropolarimetric measurements are done in only one 
single line (i.e. at one altitude), it is not possible to de¬ 
rive all the components of the electric current density but 
only the vertical one. The spatial resolution is 0.91 arcsec 
and one map of the vector magnetic field and of the verti- 
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the method lies in a better determination of the field incli¬ 
nation. This is of major importance for our study since the 
currents are determined from the horizontal components of 
the magnetic field. 

After the inversion, the 180° remaining azimuth ambi¬ 
guity was resolved by applying the ME Q code developed 
by Metcalf, Leka, Barnes and Crouch (|Leka et al.l I2QQ9I) 
and available at http://www.cora.nwra.com/AMBIG/. Af¬ 
ter resolving the ambiguity, the magnetic field vectors were 
rotated into the local reference frame (i.e. in the helio¬ 
graphic reference frame, see figure EH), where the local 
vertical axis is the Oz axis (perpendicular to the solar sur¬ 
face at the center of the frame). The vertical component of 
the electric current density (perpendicular to the plane of 
photosphere, i.e. in the heliographic coordinates) was then 
calculated via the curl of the magnetic field. 

2.2. Flare observations 



Vertical current density - 15-Feb-2011 01:48:00 



_ X (arcsec) _ 

BUI ., , _ 

-300 -200 -100 0 100 200 300 

vertical current density (mA/m^) 

Fig. 1. Top: vector magnetic field map, bottom: vertical electric 
current density map of the active region NO A A AR 11158; pro¬ 
duced from HMI data, on February 15 2011 at 01:48:00 UT, and 
represented in the plane-of-sky (see text and Appendix O for 
more details). The rectangular box represents the field-of-view 
used in the present study. Top: the arrows represent the hori¬ 
zontal component of the field (for horizontal component greater 
than 100 G), and the colors represent the intensity of the ver¬ 
tical magnetic field (see scale). Bottom: vertical component of 
the electric current density (see color scale). 

cal electric current can be calculated at a 12 minute-time 
cadence. 

The details of the calculation of magnetic f ield and verti¬ 
cal de nsity for this set of data is described in I Janvier et al.l 
(|2Q14l ). Level-lb IQUV data were used for inversion with 
th e Milne-Eddington inv ersion code UNNOEIT presented 
in iBommier et al] (|2QQ7f ). The specificity of UNNOEIT is 
that, in order to take into account the unresolved magnetic 
structures, a magnetic filling factor is introduced as a free 
parameter of the Levenberg-Marquardt algorithm that fits 
the observed set of profiles with a theoretical one. How¬ 
ever, for further application we use only the averaged field, 
i.e. the product of t he field by the magne tic filling factor, 
as recommended by iBommier et al.l (|2QQ7[ ) . The interest of 


RHESSI {Reuven Ramaty High Energy Solar Spectroscopic 
Imager) measures hard X-ray and gamma-ray emissions 
from the Sun in the 3 keV - 10 MeV range (|Lin et al.l 
|2QQ2[ ). Equipped with 9 rotating collimators (pairs of iden- 
tical grids in front of each detector), it pro vides images 
(iHnrford et al.l [200^ in addition to spectra ([Smith et al.l 
|2QQ2[ ). The spatial resolution is determined by the pitch of 
the grids of the collimators used to reconstruct the image. 
In this paper, images were obtained using the CLEAN and 
Visibility Eorward Eit algorithms using all the collimators 
except the first one (smallest pitch): the spatial resolution 
is then of 3.9 arcsec. 

The Al A Atmo spheric Imaging Assembly instrument 
([Lemen et al.ll2Ql^ aboard SDO {Solar Dynamics Observa¬ 
tory) provides images of the solar chromosphere and corona 
in EUV in six different wavelengths. The spatial resolution 
is 1.5 arcsec and for each wavelength, the time cadence is of 
12 seconds, which enables the detailed study of the evolu¬ 
tion in the flare of magnetic structures and flaring plasma. 
As the X2.2 flare on Eebruary 15 2011 is very intense, most 
of the AIA channels are saturated around the time of the 
peak of the flare, even with the minimal exposure time. 
Only the channels at 94 A and 335 A are not saturated at 
the time of the flare and therefore used for comparison with 
X-ray and magnetic observations. 

While magnetic fields and current densities are calcu¬ 
lated in the heliographic coordinates, RHESSI and AIA 
images are obtained in the plane-of-sky. Therefore, to com¬ 
bine e.g. magnetic field maps with e.g. RHESSI images, it is 
necessary to represent the magnetic fields an d cu rrent den¬ 
sities in the plane-of-sky coordinates. Eigure [AT] shows the 
geometry associated with the different measurements and 
the change of coordinates to compare the different observa¬ 
tions is described in Appendix [Aj One map of the full vector 
magnetic field and of the vertical currents represented in the 
plane-of-sky is shown in figure [T](note that the components 
of the field and vertical current density are not changed in 
this process). 

3. Observational Results 

3.1. Temporal and spatial evolution of the X-ray emissions 

The X2.2 flare on Eeburary 15 2011 was detected by 
RHESSI in the 3-100 keV energy range. The RHESSI cor- 
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Fig. 2. RHESSI corrected count rates between 01:45 and 02:15 
UT, for different energy ranges (green: 12-25 keV, cyan: 25-50 
keV, yellow: 50-100 keV) and GOES flux between 1.0 and 8.0 
A (dashed line). The vertical dashed-dotted lines at 01:48 and 
02:00 UT represent the time of the two magnetic maps used in 
the study, and the five vertical dotted lines represent the mean 
times when X-ray images are produced. 


rected count rates are presented in figure [21 together with 
the X-ray flux from GOES. In this figure the count rates 
are corrected from the changes of attenuator state during 
the flare; indeed, there was no attenuator (attenuator state 
AO) until around 01:47:10 UT, then one attenuator (Al) 
was in front of the detectors, and finally at 01:49:40 UT 
two attenuators (A3) were in place. At 02:08:50 UT, only 
the first attenuator was left in place. The peak of the X-ray 
emission in the 12-25 keV range is around 01:54 UT. 

While the temporal evolution of the count rate between 
12 and 50 keV is relatively smooth, the emission between 50 
and 100 keV is more structured. Nevertheless the different 
peaks of emission in the 50-100 keV band have counterparts 
visible in the 25-50 keV range. Images created in the 12-25 
keV, 25-50 keV and 50-100 keV energy bands, for different 
peaks of the 50-100 keV count rate, are shown in figure |3l 
We used an integration time of 16 seconds, and detectors 
2F to 9F. The first set of images (first time interval) is in 
the attenuator state Al, whereas the other images are in 
the attenuator state A3. 

Figure [3] shows that the X-ray emission sources both in 
the 12-25 keV and 25-50 keV energy ranges (in green and 
blue) have an elongated shape (more than 30 arcsec long) 
with several local maxima in space (the two main max¬ 
ima are clearly visible in the images at 25-50 keV for the 
last four time intervals in figure [3]). The 12-25 and 25-50 
keV emissions keep this elongated form during most of the 
flare. Moreover, the 25-50 keV emission is more extended 
in size than the 12-25 keV one, in each time interval. This 
has been verified first by measuring the length of the X- 
ray sources using the 50 % contours in the CLEAN images. 
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Fig. 3. Left: RHESSI images at 12-25 keV (green), with con¬ 
tours at 12-25 keV (green) and 50-100 keV (yellow) overlaid on 
the images. Right: RHESSI images at 25-50 keV (blue), with 
contours at 25-50 keV (blue) and 50-100 keV (yellow) overlaid 
on the images. Images and contours are shown for five time in¬ 
tervals integrated over 16 seconds (from top to bottom: 01:49:00- 
01:49:16, 01:51:00-01:51:16, 01:52:42-01:52:58, 01:55:02-01:55:18 
and 01:55:22-01:55:38 UT) corresponding to some peaks in the 
50-100 keV range. Images are obtained using detectors 2F, 3F, 
4F, 5F, 6F, 7F, 8F and 9F and the GLEAN algorithm. The 
green, blue and yellow contours are 60%, 70%, 80% and 90% of 
the image maximum value. 


as illustrated in figure |4] (left). To better quantify the size 
of the X-ray source, the visibility forward fitting technique 
was further used ( seelSchmahl et al']l2QQ7l for the definition 
of visibilities and IXu et al.l l2QQ8l for examples of applica¬ 
tion). This algorithm allows to compute not only the size 
of the source but also the error on the size once a model 
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Fig. 4. RHESSI images and contours at 12-25 keV (green) and 
25-50 keV (cyan) obtained between 01:50:48 and 01:51:28 UT 
using two different methods. Left: CLEAN images and contours 
at 50%, 70%, and 90% of the maximum. Right: CLEAN im¬ 
ages and contours reconstructed with the visibility forward fit 
technique (50%, 70% and 90% of the maximum). 


A more noticeable change of configuration of the 50- 
100 keV sources occurs between intervals (c) and (d) when 
sources A and B/B’ disappear and two new sources, E and 
D’ appear. The two sources furthermore seem to be foot- 
points of a loop perpendicular to the elongated structure 
seen at 12-50 keV. It can be noted than source E at 50-100 
keV is cospatial with some local maximum emission at 12-25 
and 25-50 keV and that at the same position of source D’, 
a source D also appeared in interval (b), but in the 25-50 
keV energy range. 

The spectral analysis (performed using the two func¬ 
tions vth and thick2 in the RHESSI spectral analysis soft¬ 
ware) shows that the emission is a combination of a ther¬ 
mal component with a temperature in the 25-35 x 10^ K 
range and of a non-thermal component produced by en¬ 
ergetic electrons. There is a hardening of the non-thermal 
electron spectrum starting from interval (d), as the electron 
spectrum spectral index S evolves from [-7,-6.5] in intervals 
(a), (b) and (c) to [-5,-4.6] in intervals (d) and (e). For the 
time intervals of the images shown in figure [31 the spectral 
analysis shows that while most of the emission in the 12- 
25 keV range is of thermal origin, the emission in the 25-50 
keV range is mostly produced by non-thermal electrons (the 
proportion of non-thermal emission in this energy range is 
greater than 75% in all images). 


Table 1. Loop lengths (in arcseconds) inferred from the visibil¬ 
ity forward fit algorithm for 40 second time intervals centered 
around (b) and (c). 


12-18 keV 18-26 keV 26-40 keV 40-60 keV 

(b) 28.8 ± 1.4 32.6 ± 1.0 40.4 ± 1.8 44.5 ± 4.3 

(c) 28.8 ±0.85 29.4 ± 1.3 46.3 ± 1.9 48.5 ± 3.4 


has been chosen. In order to get sufficient statistics, images 
were created with a 40 seconds integration time (instead 
of the 16 seconds of figure [3]). The result of the visibility 
forward fi t reconstruction using a curved elliptical Gaus¬ 
sian loop (|Xn et al.l 120081 ) as a model is shown in figure 0] 
(right column) for the interval 01:50:48 UT - 01:51:28 UT 
centered around interval (b). To better determine the evo¬ 
lution of the length of the X-ray source with energy, images 
were built in narrower energy bins. It is finally found that 
for the intervals centered around (b) and (c), the length of 
the source is increasing logarithmically with energy with a 
slope of 0.23 ±0.03 and 0.60 ±0.20 (see table [1]): as already 
observed in some flares (|Xn et al.ll2008[) . 

At higher energies (50-100 keV) the X-ray emission 
comes from compact sources: A and B visible in intervals 
(a), (b) and (c) and sources D and E in intervals (d) and 
(e). In intervals (a) to (c), sources A, B and B’ are localized 
along the elongated structure seen at lower energies. While 
sources B and B’ seem to be located at a footpoint of the 
elongated sources observed at 12-50 keV, source A is located 
close to the middle of this structure. Between intervals (a) 
and (c), it is noticeable that source B has moved of about 
9 arcseconds to the north-west (source B’), while source A 
has itself moved of around 3 arcseconds towards the west. 


3.2. Comparison with the EUV data from SDO/AIA 

Emissions in both 94 A and 335 A channels at around 
01:51 UT are shown in figure [5] (middle and top pannels 
respectively). While the 94 A emission shows an elongated 
structure together with more compact sources at the edge, 
at 335 A most of the emission arises from the compact 
sources. Figure [5] (middle) furthermore shows that some 
structures are emitting in both 94 A and 335 A channels 
(yellow contours), whereas other are emitting only at 94 A 
(red contours). This can be interpreted considering the dif- 
ferent channe l respo nse to source temperature. Table 1 from 
iLemen et al.l (|2Q12[ ) gives the characteristic emission tem¬ 
peratures of the plasma for the 94 A and 335 A channels: 
6.3 X 10^ K an d 2.5 x 10^ K r e specti vely. However, as shown 
in figure 13 of iLemen et al.l (|2Q12[) . the different channels 
are sensitive to a wide distribution of temperatures. From 
this figure, it is deduced that, if a structure is bright at 
94 A and 335 A, the plasma temperature is most probably 
around 1 x 10^ K (yellow contours), while if it is bright only 
in the 94 A channel, it would be much higher (more than 
6 X 10^ K, red contours). The superposition of the elongated 
part of the X-ray sources below 50 keV with the elongated 
structure seen at 94 A in figure [5] (bottom) strongly sup¬ 
ports this interpretation since the spectral X-ray analysis 
shows that the X-ray emission below 50 keV is partly emit¬ 
ted by a hot plasma with a temperature in the range 25-35 
xlO^ K. Figure [5] (bottom) also shows that source D in im¬ 
age (b) in figure [3] is roughly co-spatial with a EUV source 
detected both at 335 and 94 A. This source is most probably 
the footpoint of coronal loops in which energetic electrons 
are injected. Figure [6] shows the evolution with time of both 
X-ray and EUV sources at 335 A. Only the 335 A can be 
used for this comparison since the emission is saturated at 
94 A after 01:52 UT. In the last two images (intervals d and 
e) in which the configuration of HXR sources at 50-100 keV 
is dramatically changed compared to interval (c), it should 
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3.3. Comparison with the vertical electric current densities 

Figure 0 shows the combination of the vertical magnetic 
field and of the vertical currents for the field of view indi¬ 
cated by the black box in figure [U i.e. the part where the X- 


Fig. 6. SDO/AIA images at 335 A (grey scale), of a part of the 
active region 11158, on February 15 2011. RHESSI X-ray con¬ 
tours 12-25 keV (green), 25-50 keV (cyan) and 50-100 keV (yel¬ 
low) at 01:51:08, 01:52:10, 01:55:10 and 01:55:30 UT are overlaid 
(contour levels are the same as in figure (3]). 
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Fig. 5. SDO/AIA image at 335 A (top), and 94 A (middle and 
bottom), on February 15 2011, around 01:51:15 UT. In the mid¬ 
dle image, the region bordered in red is not clearly detected at 
335 A (most probably the plasma temperature is higher than 6 
xlO® K; see text); regions bordered in yellow are seen in both 
channels (probably the plasma temperature is close to 1 xlO® 
K). Bottom: image at 94 A on which RHESSI contours are over¬ 
laid. The green, blue and yellow contours are the RHESSI emis¬ 
sions at 12-25 keV, 25-50 keV and 50-100 keV respectively, in¬ 
tegrated between 01:51:00 and 01:51:16 UT. The contour levels 
are similar to the ones in figure [3l 


be noted that new arcade-like structures appeared in EUV 
between the two 50-100 keV sources. 
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ray emissions of the flare come from. The improved spatial 
resolution of magnetograms allows to distinguish very fin e 
structures such as the current ribbons (| Janvier et al.|[2QT3) , 
which are visible on figured] as the alignement of close posi¬ 
tive and negative currents and which are shown by the black 
arrow in the top pannel of figure [71 These current ribbons 
are present in both maps at 01:48 and 02:00 UT. Note that 
these narrow current ribbons have a fragmented appearance 
with many knots of high current densities distributed along 
the ribbons. The other noticable feature in figure [7] (bot¬ 
tom) is the appearance, in the region marked by the red 
box, of strong currents. The current density in this region 
is increased by 15% between the two maps. This evolution 
of the currents has alreadv been noted bv [Janvier et al.j 

(Eoil. 
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Fig. 7. Magnetic field maps (in grey scale) from SDO/HMI of a 
part of the active region 11158, on February 15 2011, at 01:48:00 
UT (top) and 02:00:00 UT (bottom), represented in the plane- 
of-sky. The orange and blue contours represent the positive and 
negative vertical electric current densities respectively, with the 
magnitude greater than 100 mA/m^. The black arrow in the top 
pannel indicates the current ribbons extending from 195 arcsec 
to 230 arcsec in the x-direction, and laying between -220 arcsec 
and -230 arcsec in the y-direction. The red box enlightens the 
primary difference between the two maps : in this box, the total 
negative vertical current density varies from —974.9 x 10^ A to 
-1062 X 10® A, and the total positive vertical current density 
varies from +616.9 x 10® A to +762.4 x 10® A. Therefore, the 
total increase of vertical current density is of 232 x 10® A for the 
area of the red box. This represents an increase of 15% of the 
total current density in this area. 


Figure [8] shows the evolution of X-ray sources with time, 
overlaid on the vertical magnetic field and current den¬ 
sity maps at the photospheric level obtained at respectively 
01:48:00 UT (left) and 02:00:00 UT (right). Two main re¬ 
sults can be drawn from this figure: 

— Part of the elongated X-ray sources observed between 
12 and 50 keV in intervals (a),(b),(c) overlay the current 
ribbons. 

— The new X-ray source D’ appearing at 50-100 keV in 
intervals (d) and (e) is located in the region where new 
vertical photospheric electric currents appeared between 
01:48 and 02:00 UT. 

The distance of the centroid of each X-ray source seen in 
Figure [HI to the nearest current concentration has been es¬ 
timated. Most of the HXR centroids are found to be within 
a distance lower than 4 arcseconds to the closest part of a 
current ribbon. However, two sources make exception: the 
centroid of source E (see images (d) and (e)) which is dis¬ 
tant of 8 to 12 arcsec to the nearest current ribbon, and the 
centroid of source B on image (b) which is at 6 arcsec to 
the nearest current ribbon. Figure [9] shows the X-ray images 
reconstructed with the visibility forward fit technique su¬ 
perimposed on the magnetic field and current density maps 
for intervals (b) and (c). The distance of the maximum of 
the elongated source (in the narrow energy ranges used in 
section ism see e.g. table [1]) to the nearest current concen¬ 
tration has also been estimated. This distance is found to 
decrease with photon energies, being of the order of 8 to 10 
arcseconds for energies below 26 keV and below 6 arcsec for 
energies above 26 keV. 

4. Discussion and interpretation 

4.1. How to interpret the spatial configuration of X-ray and 
EUV sources? 

During the fiare, the X-ray sources below 50 keV have an 
elongated shape (more than 30 arcseconds in length), cospa- 
tial with the EUV emission at 94 A. A large proportion of 
the X-ray emission in the 25-50 keV range is however mostly 
non-thermal (more than 75% of the total emission in this 
energy range). It can be furthermore noted that the length 
of the source is systematically increasing with energy from 
12-18 keV to 40-60 keV (see table [1]). 

Such an elongated shape of the hard X-ray emission as 
well as the increase of the source size with energy is not the 
standard situation. However, a few events have been re¬ 
ported in which non-thermal HXR emissions at 25-50 keV 
had an elongated shape resulting from thick-target inter¬ 
actio ns in a dense coronal loop (see e.g. IVeronig fc Brown! 
|2004| ) . The evolution with ener gy of hard X-ray source sizes 
was furthermore estimated bv IXu et al.l (|2008[ ) for several 
limb events on extended sources. They found, as in the 
present case, a slow increase of the source size with en¬ 
ergy (the logarithmic slope found for several events is in the 
range 0.11 ±0.04- 0.76±0.03). They also predicted the evo¬ 
lution of the source size with energy when X-ray emission 
results from a thermal population of electrons produced at 
the apex of the loop or when X-ray emission results from 
non-thermal electrons accelerated in an extended region at 
the loop apex. They found that in the case of thermal emis¬ 
sions, the source size should decrease with energy, while it 
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Fig. 8. Magnetic field maps (in grey scale) from SDO/HMI of a part of the active region 11158, on February 15 2011, at 01:48:00 
UT (left) and 02:00:00 UT (right). The orange and blue contours represent the positive and negative vertical electric current 
densities respectively, with amplitude >100 mA/m^. The green, cyan and yellow contours are the X-ray emissions (from RHESSI) 
at 12-25 keV, 25-50 keV and 50-100 keV respectively, integrated between (a) 01:49:00 and 01:49:16, (b) 01:51:00 and 01:51:16, (c) 
01:52:42 and 01:52:58, (d) 01:55:02 and 01:55:18, and (e) 01:55:22 and 01:55:38 UT. X-rays have been imaged with the algorithm 
CLEAN, using collimators 2, 3, 4, 5, 6, 7, 8 and 9. The contours corresponds to 60, 70, 80 and 90 % of the maximum for the X-ray 
emissions. 


should increase in the case of non-thermal emissions. Com¬ 
bining these predictions with the present observations of X- 
ray source sizes (systematic increase with energy from 12-18 
keV to 40-60 keV), we therefore conclude that in this flare 
the X-ray emissions in the 25-50 keV range coming from the 
elongated sources arise primarily from non-thermal emis¬ 
sions (in agreement with t he spectra l anal ysis of section 
IQ) . Furthermore following IXu et aJ] (|2008[ ) for the inter¬ 
pretation of the evolution of the source size with photon 
energy, it can be deduced that the acceleration region is 
extended: the size is found to be around 28 arcseconds and 
the density of the loop of the order of 3 x 10^^cm“^. Such 
values are in good agreement with what has been found for a 


few events \w Kontar et al. (|2Qll[ ) an d lGim et, a]J(l2ni2h . Tn 
their work. H^ et alJ ( 2008 1 modelled the transport of elec- 
trons on the basis of a standard coll i sional transport model. 
It was later shown by iKontar et al.l (|2Q14l ) that an increase 
of X-ray source size with energy can also be predicted in 
the context of a diffusive transport model of electrons, in 
which the turbulent pitch angle scattering of energetic elec¬ 
trons on magnetic fluctuations leads to an enhancement of 
the coronal HXR source relative to the footpoints. Compar¬ 
ing the predictions of this last model with our observations, 
the size of the acceleration region is found to be 22 arcsec 
(1.6 X 10^ km) and the electron mean free path of the or- 
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Fig. 9. Magnetic field map (in grey scale) from SDO/HMI of a 
part of the active region 11158, on February 15 2011, at 01:48:00 
UT. The orange and blue contours represent the positive and 
negative vertical electric current densities respectively, with am¬ 
plitude >100 mAlw?. The green, cyan, yellow and red contours 
are the X-ray emission (from RHESSI) at 12-18 keV, 18-26 keV, 
26-40 keV and 40-60 keV respectively, integrated between (b) 
01:50:48 and 01:51:28 and (c) 01:52:30 and 01:53:10 UT. X-rays 
have been imaged with the Visibility Forward Fit technique us¬ 
ing collimators 2, 3, 4, 5, 6, 7, 8 and 9. The contours corresponds 
to 50, 70 and 90 % of the maximum for the X-ray emissions. The 
crosses are at the location of the maxima of the sources. 


der of 10^ c m, whi ch is in agreement with the results of 
iKontar et al.l (|2Q14l ) for another flare. 

At higher energies (>50 keV), the X-ray emission is pro¬ 
duced by a non-thermal population of electrons and the 
emissions come from more compact sources: sources A, B 
and B’ in images (a) to (c) and sources D’ and E in images 
(d) and (e). While sources B and B’ appear to trace the 
impact of energetic electrons at the footpoint structures, 
the 50-100 keV source labelled A in images (a), (b) and (c) 
is located in the middle of the elongated structure observed 
at lower energies. Its location is also very close to the posi¬ 
tion of the maximum of the elongated X-ray source in the 
40-60 keV range. This high-energy source could be similar 
to coro nal sources more eas i ly observed in the case of limb 
events (Kmcker et al.ll20()^. l2010l : IChen fc Petrosiai]l2012l: 
ISn et al.ll2013[ ). but projected on the disk in our case. The 
combined location of coronal and loop-top sources could 
trace the location where the main reconnection process oc¬ 
curs i.e. above the loop in which electrons then propagate 
(see e.g. ISn et al.ll2Q13[) . In this context, the evolution of 
the location of sources A, B and B’ between time inter¬ 
vals (a) and (c) most probably follow, the evolution of 
the acceleration site. It must be recalled that movement 
of X-ray footpoints are usually observed in the course of 
flares and related to the magnetic reconnection process (see 


e.g. iFletcher et al.l l20lll for a review). A more significant 
change in the configuration of HXR sources above 50 keV 
is observed later in the flare (images (d) and (e)). Simul¬ 
taneously with the hardening of the non-thermal electron 
spectra, new high energy HXR sources appear (sources D’ 
and E) most probably located at the footpoints of an ar¬ 
cade of loops which appeared in EUV at the same time. 
These observations show that after « 01h55 UT, magnetic 
reconnection and the resulting plasma heating and parti¬ 
cle acceleration occurs in a new structure. It is finally no¬ 
ticeable that the position and orientation of the new EUV 
bright structure are consistent with the ones of the post¬ 
flare loops ap pearing in th e 3DMHD simulation of the flare 
performed by llnone et al.l (|2Q14l ) . 

4.2. How to interpret the relation between X-ray emissions 
and photospheric vertical electric currents? 

As discussed in section loi part of the elongated thermal 
and non-thermal X-ray sources (observed at energies be¬ 
low 50 keV) are located above the narrow and elongated 
current ribbons measured at the photospheric layer. These 
elongated X-ray structures have furthermore been identi¬ 
fied in the previous section as being produced in the corona 
by electrons injected from an extended acceleration region 
located close to the apex of the magnetic structures. At 
energies above 50 keV, the X-ray source A is furthermore 
located very close (in projection) to the main current rib¬ 
bon (shown by a black arrow in the top panel of Figure Q 
revealing the close association between the location of the 
acceleration region and the location of the electric current 
sheets. 

As it has been discussed in I Janvier et al.l (|2014l ). the 
photospheric current ribbons (closely related to the flare 
EUV ribbons) are the tracers at the photospheric boundary 
of the electric current sheets present in the coronal volume. 
These current layers are themselves formed in regions of 
strong gradients of ma gnetic connectivit y (see figure 7 of 
iJanvier et al.l l2Q14l and IZhao et al.l l2Q14l for the computa¬ 
tion of QSLs in the region) known to be th e preferred loca¬ 
tions where reconnection can occur (see e.g. lDemonlin et aJ] 
Il996[ ). As a consequence, plasma heating and particle ac¬ 
celeration can be produced in these regions. 

This explains why part of the thermal and non-thermal 
X-ray emissions produced in the coronal elongated sources 
are found above the photospheric current ribbons which 
trace the footprints of the coronal current layers. 

At energies above 50 keV, the apparition of the X-ray 
source D’ is accompanied with the appearance of cospa- 
tial photospheric currents in the same time interval. Note 
that the observation of the current density evolution has 
to be carefully examined since polarization measurements 
during flares can be biaised, and some of the observed fea¬ 
tures could result from polarization artifacts due to im¬ 
pact of non-thermal particle beams (see lHenonx fc Karlick\^ 
l2Q13l. and references therein) or induced by resonant scat¬ 
tering d ue to radiation anisotropies at the edges of flare 
ribbons ([Stepan &: Heinzelll2Q 131) . [Janvier et al.l (|2Q14l ) dis¬ 
cussed this issue for this specific event. They analysed cau¬ 
tiously the vector magnetic field in the active region and 
concluded that the changes in the electric current density 
are not artifacts due to polarization effects. They argue 
that the horizontal fields (from which vertical electric cur¬ 
rents are derived) are consistent in both maps at 01:36 
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01:48 UT 


02:00 UT 


02:48 UT 



Fig. 10. Vertical current density (color) and horizontal magnetic field (arrows, see scale), for the region indicated by the black box 
in figure [7| (28" x 26" centered at x=195 arcsec and y=-214 arcsec), at 01:48 UT (left), 02:00 UT (middle) and 02:48 UT (right). 
The region where the current density increases by 15% between 01:48 and 02:00 UT is surrounded by the black rectangle on the 
three present figures. 


UT and 02:00 UT with a signal free of artifacts given the 
smooth rotation of the magnet ic field in this regio n. Fig¬ 
ure [10] (inspired from figure 4 in I Janvier et al.|[2014[ ) clearly 
shows a coherent direction of the horizontal magnetic field 
which varies smoothly within the field of view. Moreover, 
the curvature of the horizontal field lines increases between 
01:48 and 02:00 UT, which is responsible for the increase 
of the observed vertical current (see also the observations 
bv lPetrijl2013[ ). Finally, the increase in the current density 
shown in figure 018 persistent, as it can be seen in figure [TUI 
This is a very strong argument against artifacts due to en¬ 
ergetic particle precipitation or flare radiation anisotropies, 
because at 02:48 UT the GOES X-ray flare has largely de¬ 
cayed . We thus reach the same conclusion as I Janvier et al.l 
(|2014[ ): the increase of the photospheric currents in the 
black box of figure [10] is not due to a polarization artifact 
but reflects real changes in the horizontal magnetic field 
and is thus real. 


As a conclusion, the appearance of the new HXR source 
D’ and the increase in the same region and in the same 
time interval of the photospheric currents are linked. The 
increase in the photospheric currents can be interpreted as 
a response at the photospheric layer of the change of mag¬ 
netic topology and current systems in t he corona induced 
by reconnection (see 1 Janvier et alll2Q14l) . As particle accel¬ 
eration is also a consequence of the magnetic reconnection 
in the corona, it naturally explains the close temporal and 
spatial association between locally enhanced X-ray emis¬ 
sion and increase in the same region and at the same time 
of photospheric currents: both result from the same phe¬ 
nomenon, namely magnetic reconnection in the localized 
coronal current sheets. It should finally be noted that such 
a related evolution of X-ray sources and electric currents 
in the course of a flare has not been reported previously, 
given the low cadence of magnetic field measurements be¬ 
fore SDO/HMI. However, even in the present case, it is 
clear that the cadence to derive magnetic field maps (and 
current density maps) from HMI (12 minutes) is still too 
low to be able to track the evolution of electric currents on 
the time scales relevant to the evolution of X-ray sources 
(a few tens of seconds). 


4.3. Comparison with previous results 

As recalled in the introduction, the link between vertical 
electric currents and electron precipitation sites as diag¬ 
nosed by HXR emissions at footpoints has been investigated 
in the 1990s combining YOHKOH/HXT observations and 
vecto r magnetograph data from the M ees Solar Observa¬ 
tory (|Li et al.lll997l : [infield et al.lll992f ). It was found that 
the electron precipitation sites as detected via the HXR 
footpoints are not exactly co-spatial with the regions of 
highest vertical current densities but adjacent to the cur¬ 
rent channels. In this new study, using vector magnetograph 
data (and thus vertical electric maps) obtained from space 
(with an improved seeing) with a higher spatial resolution 
in combination with HXR images obtained at the same time 
and with improved image dynamics and spatial resolution, 
the link between electron precipitation sites and current 
concentrations has been revisited. These new observations 
confirm some of the previous result that HXR footpoints 
indicating electron precipitations sites are not exactly cos¬ 
pat ial with the region of highest electric current densities. 
This is in particular the case in this flare of the footpoint 
sources B and E which are at a distance of 6 to 12 arcsec 
to the nearest current ribbons in the present flare. 

HXR footpoints can also be found at the same loca¬ 
tion as the current ribbons themselves. This is the case of 
the HXR footpoint D’ seen at energies above 50 keV af¬ 
ter 01:55 UT which is cospatial with new current ribbons 
which appeared themselves between 01:48 and 02:00 UT. It 
should be emphasized that the observation of such a com¬ 
bined evolution of HXR footpoints and electric currents is 
reported here for the first time. A complete interpretation 
of this observation is not within the scope of this paper 
aimed at observations but would require a more complete 
study combining models of 3D reconnection in flares and 
particle acceleration. 

In addition to HXR footpoint sources, a significant part 
of the HXR emissions observed in this flare arises from an 
extended coronal source observed up to 60 keV. Such coro¬ 
nal HXR sources were not observed in the events studied 
by (|Li et al.lll997l:[infield et al.lll992[ ). The comparison of 
coronal HXR sources with electric currents ribbons is then 
performed for the first time in this study. It is found that 
part of the thermal and non-thermal coronal sources just 
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overlay the current ribbons. This should be further inves¬ 
tigated in other events and interpreted in the context of 
combined models of 3D reconnection in flares and particle 
acceleration but this clearly provides support to the mod¬ 
els of p article acceleration in reconnecting current sheets 
(see e.g. IZharkova et al.|[2Qlll. for a review of acceleration 
models). 

5. Summary and conclusion 

This study is the first detailed comparison between the spa¬ 
tial distribution of photospheric vertical currents and of X- 
ray emission sites produced by hot plasma and accelerated 
electrons. For the first time, X-ray images and maps of elec¬ 
tric currents (derived from vector magnetograms) were ob¬ 
tained at the same time with an improved time cadence for 
the observations of vector magnetograms (and derivation of 
electric current maps) which furthermore enables to study 
the coupled evolution of the photospheric electric currents 
and the HXR emission sites in the course of the flare. 

The main results obtained from this first study are the 
following ones: 

— Parts of the thermal and non-thermal X-ray emissions 
produced in an extended coronal source overlay the elon¬ 
gated narrow current ribbons observed at the photo¬ 
spheric level. 

— A new hard X-ray source at energies above 50 keV ap¬ 
pears in the course of the flare in association with an 
increase of 15% of the photospheric current at around 
the same time and the same location. This shows here a 
clear link between particle acceleration and reconnect¬ 
ing current sheets. 

These two results can be qualitatively explained in the 
context of the commonly admitted scenario in which mag¬ 
netic reconnection occurs at current-carrying QSLs in the 
corona with part of the energy released transferred to 
plasma heating and particle acceleration. Since X-ray (and 
EUV) emissions are signatures of plasma heating and of 
particle acceleration, and since photospheric currents can 
trace the footprints of the coronal currents embedded in 
QLSs, some spatial and temporal correlation can be natu¬ 
rally expect between X-ray (and EUV) emitting sites and 
photospheric current ribbons. In addition, the evolution of 
magnetic reconnection sites in the course of the flare, may 
lead to two linked consequences: on one hand, plasma heat¬ 
ing and particle acceleration are produced at different lo¬ 
cations, potentially leading to the appearance of not only 
new EUV post-flare loops but also of new X-ray sources; on 
the other hand, the change of magnetic topology due to the 
evolving reconnection process leads to an increase of the 
photosp heric current densiti es as the same place and time 
(see also I Janvier et al.ll2Q14l ). These new results have been 
obtained so far on only one flare. A similar study should be 
further extended in the future to other X-class flares, to ex¬ 
amine whether this relation between photospheric current 
ribbons and X-ray sources is systematically observed.This 
work furthermore shows that the evolution of horizontal 
magnetic Adds and of derived photospheric vertical elec¬ 
tric currents can be observed on timescales of 12 minutes 
during a strong flare. The cadence is however still too low 
to derive the detailed timing between the evolution of the 
fields and currents and the evolution of the energetic elec¬ 
tron acceleration and interaction sites. Einally, this kind of 


studies would deserve future interpretation and modeling: 
in particular tentatively coupling acceleration processes to 
3D MHD models should allow to derive the location of the 
reconnecting current sheets and of their observed photo¬ 
spheric traces together with the location of HXR sources. 
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Appendix A: Spatial superposition of magnetic 
field, electric current density and emissions 
from flares 

Appendix A.1: Change of Frame 

To compare the X-ray emissions with magnetic held and 
current density maps, the frame of reference must be care¬ 
fully chosen. Indeed, the X-ray images are created in the 
cartesian coordinates (Plane-Of-Sky), whereas the mag¬ 
netic maps were taken and derivated in the heliographic 
system (heliographic plane, see figure [AT]) . We chose to use 
in the paper the cartesian coordinates, since X-ray sources 
vary in timescales of the order of seconds, whereas the mag¬ 
netic held maps are obtained at a cadence of 12 minutes. 

The relation between (x,^), the cartesian coordinates of 
a point and its heliographic coordinates (T,6) is given by: 

X = sin(L - Lc) cos(5) . . 

y = sin(6) cos(Bc) — cos(L — Lc) cos(6) sm{Bc) 

with L the longitude of the point relative to the central 
meridian, b the latitude of the point. Be and Lc the lat¬ 
itude and longitude of the center of the solar disk at the 
time of observation. 

The change of frame also modifies the pixel size. In the 
case of an active region located near the disk center, several 
approximations can be made and the pixel size is described 
by: 

pXcar - PXhelio COs(L - Lc) 

Wear - pXhelio{cOs{b) COs{Bc) (A.2) 

+ cos(L — Lc) sin(6) sm{Bc)) 

with {pXcamPUcar) the pixel size in the cartesian coordi¬ 
nates (in the POS), and {pxheiioiPVheiio) the pixel size in 
the heliographic coordinates (in the heliographic plane). 

Appendix A.2: Time gap between different observations 

The second problem when comparing X-ray emissions and 
magnetic held maps is the time gap between the different 
observations. We use the time of the X-ray images as a 
reference since X-ray sources vary faster than the magnetic 
held, and also because we generally have more images in 
X-rays than maps of the magnetic held. 

The time gap between the X-ray images and the mag¬ 
netic maps must be taken into account because the rotation 
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Fig. A.l. The heliographic coordinates: line emission used for 
spectropolarimetry comes from one layer of the Sun (circle). The 
magnetic field and electric current density are then calculated 
in the heliographic coordinates, as the drawing shows: the hori¬ 
zontal component is in the heliographic plane, which is tangent 
to the considered layer of the Sun, and the vertical component is 
perpendicular to that, along the (Oz) axis. The X-ray emissions 
are taken in the plane-of-sky (POS), i.e. in the plane perpendic¬ 
ular to the line-of-sight (LOS). 


of the solar surface is not negligible. It has been corrected 
by setting the observation time of the X-ray image when 
looking for Be and Lc in the ephemerids to use the equa¬ 
tion rvn Note that we used the get_sun procedure in the 
SolarSoft to get the ephemerids {Be and Le at time of 
observation). 
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